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The  Gurney  method  was  originally  developed  to 
estimate  terminal  velocities  Imparted  to  metals  by  ex¬ 
plosives,  and  was  subsequently  uS3d  to  correlate  ex¬ 
ploding  foil  system  performance  with  firing  circuit 
electrical  parameters.  The  Inclusion  of  Internal  en¬ 
ergy  In  the  explosive  gas  permits  calculation  of  the 
acceleration  history  of  the  system.  Comparison  with 
experimental  and  hydrocode  results  Indicates  that  pro¬ 
per  selection  of  Initial  displacement,  velocity,  and 
the  polytropic  gas  exponent  can  bring  model  and  ex¬ 
perimental  results  Into  close  agreement. 


I.  INTRODUCTION 

We  want  to  make  some  remarks  today  on  the  use  of  the  Gurney  method  to  esti¬ 
mate  the  acceleration  history  of  an  exploding-foll -driven  flyer.  Schmidt,  Seitz, 
and  Wackerle^  formulated  the  equation  of  motion  for  such  a  system  with  time  de¬ 
pendent  energy  deposition.  As  McGlaun  has  shown,  the  approximation  of  uniform 
pressure  and  density  throughout  the  foil  gas  can  be  replaced  by  a  continuum  so¬ 
lution  which  Improves  the  overall  accuracy.  This  modification  could  be  added  to 
the  work  I  am  presenting. 

Today  we  want  to  discuss  the  governing  differential  equation  (D.E.)  of 
motion  and  mention  some  of  the  things  we  can  learn  from  It.  We  want  to  discuss 
some  things  that  the  hydrocodes  can  tell  us  about  solving  the  D.E.  And,  finally, 
we  want  to  discuss  some  of  the  things  we  can  learn  from  It.  We  want  to  discuss 
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some  of  the  things  we  can  do  with  the  D.E.  that  aren't  so  easy  with  the  hydro 
codes. 


II.  THE  DIFFERENTIAL  EQUATION 

The  D.  E.  of  motion  can  be  written  In  the  form  shown  In  Figure  1,  which 
gives  an  expression  for  the  acceleration  of  the  flyer,  x.  We  see  that  the  driv¬ 
ing  term  Is  the  quantity  In  brackets.  It  Is  the  difference  between  the  Input 
energy  and  the  system  kinetic  energy  and  thus  Is  the  Internal  energy.  In  this 
loss-less  model,  we  see  all  the  Input  energy  will  be  converted  to  kinetic  energy 
In  a  long  barrel. 

The  Internal  energy  Is  multiplied  by  several  factors.  Generally  we  want 
rapid  acceleration,  so  we  want  the  multipliers  to  be  large.  The  factor  j[  Is 

large  when  flyer  displacements  are  small.  Thus,  Input  energy  Is  more  rapidly 

converted  to  kinetic  energy  If  Inputed  when  X  Is  small.  The  factor  ^  Is  the 

"charge  to  mass  ratio".  It  Is  comforting  to  see  that  acceleration  Is  greater 
for  large  charge/mass  ratios.  Likewise,  (HA)~^  Is  largest  when  A  Is  small 
corresponding  to  heavily  tamped  systems.  Finally,  the  polytropic  gas  exponent, 
n.  Is  a  factor. 

III.  HYDROCODE 

When  Schmidt,  Seitz,  and  Wackerle  carried  out  the  solution  they  used  for  n 
a  value  of  5/3  corresponding  to  the  adiabatic  gas  constant  for  a  monatomic  gas. 
This  leads  to  velocity-time  results  much  higher  than  experiment.  To  bring  re¬ 
sults  Into  agreement,  they  obtained  an  empirical  correction  factor  for  the  power 
pulse  Illustrated  schematically  In  Figure  2.  The  Implication  of  this  result  Is 
that  early  Input  energy  does  not  contribute  to  flyer  motion  while  late  time 
energy  Input  Is  used  at  full  value.  As  we  shall  see  there  are  other  means  to 
obtain  agreement  between  experiment  and  the  model. 

We  want  to  propose  an  alternate  method  for  bringing  this  model  Into  agree¬ 
ment,  and  It  Involves  selection  of  different  values  for  n  and  for  the  Initial 
conditions  for  the  D.  E. 

From  CHART  D  one-dimensional  hydrocode  simulations,  we  looked  for  pressure- 
volume  relationship  In  the  expanding  foil  gas.  CHART  D  has  carefully  constructed 
analytic  equations  of  state  for  aluminum  and  copper  which  are  valid  over  a  large 
range  of  tanperatures  and  densities  and  Include  solid,  liquid,  vapor,  mixed 
phases  and  plasma.  In  CHART  D  simulations  of  foil /flyer  systems,  the  P-V  re¬ 
lationship  after  shut-off  of  Input  energy  gave  n  •  1.23  for  alumlnun.  To  check 
this  value,  we  used  the  E.O.S.  portion  of  CHART  D  for  aluminum  which  gave 
1.3  <  n  <  1.4  In  the  range  of  temperatures  and  densities  of  Interest,  which  Is 


close  to  the  value  obtained  In  the  simulation. 
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IV.  RESULTS 


Figure  3  shows  the  effect  of  n  on  the  gas  dynamic  Gurney  velocity- sime 
calculation  for  a  flyer  using  the  power  pulse  without  empirical  correction. 

Shown  for  comparison  is  CHART  D  result. 

In  addition  to  the  polytropic  gas  constant,  the  hydrocode  furnishes  other 
useful  information  for  solving  the  D.E.  The  above  solutions  assume  Xq  =  Xo  =  0. 
However,  the  D.E.  applies  to  a  gas  dynamics  problem.  So  it  is  appropriate  to 
ask  what  are  the  initial  conditions  for  the  gas-dynamics  problem.  From  CHART  f) 
runs,  the  foil  has  been  heated  to  vaporization  at  a  time  i/hen  it  has  expanded 
about  three-fold  and  the  flyer  has  attained  about  ten  percent  of  its  impact  ve¬ 
locity.  Using  these  initial  values,  the  D.E.  solution  is  in  close  agreement 
with  CHART  D  as  shown  in  Figure  4. 

These  results  apply  to  a  small  aluminum  foil  system.  For  contrast  we  can 
also  look  at  a  large  copper  system  (Figure  5)  and  see  that  a  similar  n  is  ap¬ 
propriate  and  in  Figure  6  a  fair  agreement  can  be  produced  with  proper  initial 
conditions. 

All  of  these  results  are  aimed  at  comparison  between  two  models,  gas  dyna¬ 
mic  Gurney  D.E.  and  CHART  D  hydrocode.  As  has  been  pointed  out  before,  using 
CHART  D  with  measured  power  profiles  over-predi cts  V{t).  So,  the  D.E.  method 
would  also  overpredict.  Ue  believe  the  problem  lies  in  power  measurements  and 
we  have  an  experimental  program  underway  to  attempt  to  improve  the  measurements. 


V.  GURNEY  vs.  HYDROCODE 


In  conclusion,  let  us  comment  on  the  use  of  the  Gurney  D.E.  vs.  hydrocodes. 

We  think  the  D.E.  offers  several  useful  features. 

1.  The  D.E.  offers  Insight  into  trends  and  limiting  cases. 

2.  The  agreement  obtained  indicates  wave-propagation  effects  are  small,  so 

the  use  of  a  hydrocode  may  be  unnecessary.  , 

3.  Big  codes  on  big  computers  are  scxnewhat  intimidating.  The  D.E.  can  be 
handled  with  less  computer  expertise  and  can  be  approximated  with  tabu¬ 
lar  or  hand  calculator  techniques. 

4.  The  D.E.  is  convenient  for  use  with  fire  set  circuit  codes. 
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GAS  DYNAMICS  GURNEY  MODEL  EQUATION  OF  MOTION 
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